Thermodynamic potential of interacting Bose-Einstein gas 
confined in harmonic potential 
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We investigate the interaction effect between atoms and the finite size effect of a Bose-Einstein 
gas at finite temperature. Using a mean field approach, we derive the thermodynamic potential on 
finite systems and obtain the condensate fraction, the chemical potential and the internal energy 
with the interaction and the finite size effects taken into account. In this approach, we can treat 
thermodynamic quantities as those of the ideal Bose gas. ft is found that the interaction effect shifts 
the transition temperature by a few percent toward the lower temperature. Above the transition 
temperature, these quantities agree with the noninteracting Bose gas. 
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I. INTRODUCTION 

Bose-Einstein condensation (BEC) was first observed 
in dilute ultracold alkali atoms of rubidium Q , lithium Q 
and sodium ||. Furthermore, Mewes et al. and Enshcr 
et al. measured the condensate fraction and the energy of 
rubidium atoms . It revealed that the transition tem- 
perature To is shifted not only by the interaction effect 
but also by the finite size effect by a few percent, where T 
denotes the transition temperature of the noninteracting 
Bose gas within the external field in the thermodynamic 
limit. This means that both the interaction and the fi- 
nite size effects play an important role in the real Bose 
gases. In fact, the number of trapped atoms is typically 
N = 10 3 ~ 10 6 , and this may not be sufficiently large 
to take the thermodynamic limit. For the noninteracting 
Bose gas, the condensate fraction confined in a spherical 
harmonic potential can be analytically given as 
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where the first two terms in the right hand side corre- 
spond to the condensate fraction in the thermodynamic 
limit and the third term is the finite size correction. 

On the other hand, these experimental findings have 
stimulated much interest in the theory of the interacting 
Bose gas. Since BEC occurs when atoms are dilute and 
cold, we can treat the interaction of atoms as the two- 
body interacting Bose gas. In this case, the s-wave scat- 
tering length characterizes the strength of the two-body 
interaction. Under these conditions, the Gross-Pitaevskii 
(GP) equation can well describe the behavior of the in- 
teracting Bose gas at zero temperature [|]||. This is a 
mean field approach for the order parameter associated 
with the condensate. Using the GP equation, several au- 
thors jlO) studied the ground state and the excitation 
properties of the condensate. 

To study BEC at finite temperature, the Bogoliubov 
approximation at zero temperature was extended by Grif- 
fin El. It is called the Hartree-Fock-Bogoliubov (HFB) 



theory. In particular, Popov approximation to the HFB 
theory has been employed to explain the experimental re- 
sults [^2| . This approximation neglect the anomalous av- 
erage of the fluctuating field operator in the HFB theory. 
With the Popov approximation, thermodynamic quan- 
tities, such as the condensate fraction and the internal 
energy, agree with the experimental results at low tem- 
perature. However, recent experimental data of the exci- 
tation frequencies of the condensates indicate a discrep- 
ancy with these approximate results near the transition 
temperature T (T/T > 0.6) Q. 

In this paper, we derive the thermodynamic potential 
for the interacting Bose gas within a mean field approx- 
imation which is different from Bogoliubov approxima- 
tion. With this thermodynamic potential, we obtain the 
thermodynamic quantities; the condensate fraction, the 
entropy and the internal energy. These quantities have 
the same form as that of the ideal Bose gas, but the en- 
ergy E{ contains the interaction effect and varies as a 
function of the temperature T and the total number N. 
Therefore, we must determine not only the chemical po- 
tential fi but also the energy E{ as a function of T and 
N in a self-consistent fashion. In particular, the chem- 
ical potential determines the transition temperature on 
the condition that the chemical potential [i is equal to 
the ground state energy E Q at T = T . Here, we choose 
the repulsive interaction where the scattering length a 
is positive, and numerically solve the chemical potential 
and the energy E\. By determining the chemical poten- 
tial near T ~ T , we find that the interaction effect shifts 
To by a few percent toward the lower temperature which 
should be compared with the previous estimations of the 
semiclassical approximation ]l4| ]. 

This paper is organized in the following way. In the 
next section, we derive the thermodynamic potential for 
the interacting Bose gas at finite temperature. Then, we 
obtain thermodynamic quantities in a self-consistent way. 
In section 3, we carry out the numerical evaluation of the 
chemical potential, the condensate fraction and the inter- 
nal energy, in particular near Tq. Section 4 summarizes 
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what we have clarified in this paper. 



II. THERMODYNAMICS FOR INTERACTING 
BOSE GAS 

We derive the thermodynamic potential for the inter- 
acting Bose gas within a mean field approach. To be- 
gin with, we construct the effective Hamiltonian in this 
system. The Hamiltonian for the interacting Bose gas 
confined in an external potential V ext can be written as 



where H n (^) is the Hermite polynomial, and N n is the 
normalization constant which can be written as 
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where aho denotes the harmonic oscillator length 
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This is also an important length scale in this system. It is 
typically of the order of fim. From Eqs. ( [P| ) and (|2.6|) , 
the Hamiltonian can be written as 



±J drdr'$\r)&(r>)V(r-r>)*(r')$(r), (2.1) H = J2 £ » a » a * 



where ^(r) is the boson field operator and V(r—r') is the 
two-body atomic potential. Here, we choose the spheri- 
cal harmonic potential for the V ex t for simplicity. Thus, 
the Vext is given as 



Kxt(r) = ^muj 2 r 2 . 



(2.2) 



Since we are interested in the dilute and cold Bose gas, 
we replace V(r — r') by an effective interaction of the 
^-function type, 



K(r-r') =^ (3) (r-r'), 



(2.3) 



where the coupling constant g is related to the s-wave 
scattering length a through 
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The scattering length is an important length scale in this 
system. A typical value of the scattering length a is of 
the order of nm. The boson field operator for the ideal 
system can be written as 



*(r) 



a k e 



(2.5) 



But for the general case, the corresponding field operator 
is a sum over all normal modes 11 511 



=I>nXn(r), 



(2.6) 



when x n (r)'s define any complete set of normalized 
single-particle wave function a n is a bosonic annihila- 
tion operator for the single-particle state n. In our case, 
X n (r) is given as 
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Now we consider the first order perturbation theory 
for the thermodynamic potential (Tq |. The first order 
correction is given by the diagonal part of the interact- 
ing Hamiltonian in Eq. ( p. 10 ). This corresponds to ig- 
noring the anomalous average. Therefore, the effective 
Hamiltonian can be written as 
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Now we rewrite Eq. (2.13) in the following way 
Heff = Hq + H , 



with 
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Here, we minimize H' by taking v\ the mean occupation 
numbers 



V\ = (fli) = 
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Then, the thermodynamic potential Q = —pV can be 
given as 



n 
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where the first term in the right han d sid e is the first 
order correction. We note that Eq. ( 2.1 9| ) can be also 
determined by the condition 
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We obtain the thermodynamic qu antiti es fro m Eq . 
( |2.20D by solving Ei and fi from Eqs. fl2.18| ) and ( ]2.20| ) 
self-consistently. This is a mean field approach within 
the first order perturbation theory. The chemical po- 
tential cannot exceed the energy E\ at any temperature. 
At the transition temperature, the chemical potential is 
equal to the ground state energy. Therefore, the transi- 
tion temperature has to be determined by the che mica l 
potential. We note that the second term of Eq. (2. IS) 



shifts the transition temperature. If we choose tAj as a 
constant i.e. U, then it only causes the constant shift for 
all E{, but the transition tempera ture does not change. 
Therefore, the second term of Eq. ( [2.18 ) plays an impor- 
tant role to understand the interacti ng B ose gas. We also 
note that, when we put g = 0, Eq. (2.22) reproduces the 
total number of the noninteracting Bose gas confined in 
the spherical harmonic potential, because there are no 
constraints in between the chemical potential fi and the 
energy E\. Therefore, when we consider the interaction 
effect, we can also investigate the finite size effect at the 
same time. 



This condition means that the thermodynamic potential 
is a minimum with respect to any change of state at 

constant T, V and /i in a state of thermal equilibrium. 
Now we derive the thermodynamic quantities from 

Eqs. (2.20) a nd ( 2.21). First, we obtain the to tal nu mber 

N from Eq. fl2.20|) with the condition of Eq. fl2.2l|) 
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We can determine t he ch emical poten tial fi and the en- 
ergy Ei from Eqs. ( |2.18j ) and ( |2.22[ ) which should be 
solved self-consistently. Then, we obtain the thermody- 
namic quantities through Eq. (2.2C) with Eq. (2.21). For 



example, the entropy S can be given by 
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where we also use the condition Eq. (2.21). Finally, the 
internal energy E can be given as 



III. NUMERICAL RESULTS 



We present the numerical re sults of the rmod ynamic 
quantities which appear in Eqs. ( [2.22] ) and ( [2.24 ). First, 
we numerically solve the chemical potential [i and the 
energy E\ by using Eqs. (2.18) and ( 2.22| ) with the iter- 
ation method. Then, we obtain fj, and Ei as a function 
of T and N. If the temperature of this system is low- 
ered at fixed N, the chemical potential [i given by Eq. 
( [2^22 ) increases. It finally reaches the ground state en- 



ergy Eq, since the chemical potential cannot exceed the 
ground state energy. Therefore, we can define the tran- 
sition temperature T c as fi = Eq at T = T c . We note 
that the transition temperature T c for the noninteract- 
ing Bose gas confined in the spherical harmonic potential 
can be given as 



T c = T [ 1 
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where the second term is the finite size effect, and To 
denote the transition temperature in the thermodynamic 
limit which can be written as 



To = — 
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Because the interaction effect also shifts the transition 
temperature T c , we have to redefine the transition tem- 
perature. 
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For T > T , we solve Eqs. fl2.18| ) and ( |2.22| ) with 
fixed iV and given T. Here, the total number N and 
the strength of the interaction with the repulsive force 
are chosen to be N = 1000 and a/a/ lo = 7 x 10~ 3 , re- 
spectively. In Fig. [I], we show the calculated values of 
the chemical potential /i/fe^To as a function of T/Tq. As 
can be seen, when T/Tq ~ 0.91, /i//i;_bTo approaches to 
a finite value. This means that the chemical potential (i 
reaches the ground state energy Eq. Therefore, the tran- 
sition temperature becomes T/Tq ~ 0.91 (see, Fig. ||). 
Comparing this result with the transition temperature 
T c /T ~ 0.93 given by Eq. flo] ), the interaction effect 
further shifts the transition temperature about 2 percent 
toward the lower temperature. 

For T < Tq, the chemical potential fj, is equal to the 
ground state energy Eq. In this case, since the ground 
state occu pation numbers Nq diverges, we have to rewrite 
Eq. Q2.18D as follows 



Ei = Si + Two U ifi N 

O-ho 



aho 
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iVo+Sj^o( n j)- Under this 



where we use the relation N 

condition, we solve Eq. flOfl with n/k B T a = 0.178. We 
note that h/UbTq is 0.159 for the noninteracting model 
at the transition temperature, since /i is 3/27iu; in our 
approach. 

In Fig. ||, we show the condensation fraction as a func- 
tion of T/Tq. The dotted line refers to the condensation 
fraction of the noninteracting Bose gas in the thermo- 
dynamic limit, and the solid line represents the conden- 
sation fraction including the finite size effect (see, Eq. 
(LI)). Finally, the condensation fraction with the inter- 
action and the finite size effects is shown as solid circles. 
These correspond to each other at T > To. This means 
that we can ignore the atom-atom interaction because 
the number density is low at T > To. On the other hand, 
at near the transition temperature, atoms must occupy 
the ground state energy Eq. Since the number density 
becomes high at this temperature, the atom-atom inter- 
action becomes strong, and shifts the transition temper- 
ature. Although the finite size effect can also shift the 
transition temperature, it disappears at large N limit. 
Therefore, the interaction effect is more important than 
the finite size effect. 

In Fig. [|, we show the results for the internal energy 
E/NksTo as a function of T/Tq. It reveals that the inter- 
nal energy curve changes slope near the transition tem- 
perature. We note that the specific heat can be calcu- 
lated by differentiating the internal energy with respect 
to the temperature T. Therefore, this result means that 
the specific heat is discontinuous at the transition tem- 
perature. As can be seen from the result of Ensher et al. 
H , the derivative of the internal energy is discontinuous 
at the transition temperature. This behavior of our cal- 



culated internal energy is in a good agreement with the 
experimental result. 

Finally, we compare our results with the semiclassical 
results. In the semiclassical approximation, the dimen- 
sionless parameter Ei n t/Eki n = Na/aho is important to 
discuss the two-body atomic interaction energy E int com- 
pared to the kinetic energy Eki n ■ For the Thomas-Fermi 
(TF) limit (Na/aho S> 1), the chemical potential hq can 
be written in terms of this parameter ]l7j , |T8"| 



Hlu ( a 
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(3.4) 



In our case, Ho/ksTQ is 0.342. This is about twice as 
large as the /x//cbT = 0.178. Since we take the pa- 
rameter Na/aho <~ 1 , it does not achieve the TF limit. 
Therefore, the kinetic term plays an important r ole in 
this value. Further, the interaction term Uij in Eq. ( 2. IS ) 
suppresses the total number N of Na/aho- The shift of 
the transition temperature STq/Tq in the semiclassical 
approximation for the large TV limit can be given as um 



To aho 



(3.5) 



In our case, the shift of the transition temperature STq/Tq 
is about 3 percent. Because of the interaction term 
the transition temperature is also suppressed by its term 
in our model. 



IV. CONCLUSIONS 

We have presented the thermodynamic potential of the 
interacting Bose gas confined in the spherical harmonic 
potential within a mean field approach. In this approach, 
we minimize the interacting Hamiltonian H' under the 
condition that the mean occupation numbers v\ should 
be the Bose distribution function. Then, wc can derive 
the thermodynamic potential in a self-consistent fashion. 
With this thermodynamic potential, we obtain the ther- 
modynamic quantities which have the same form as that 
of the ideal Bose gas. Further, we carry out the numeri- 
cal calculations of the chemical potential, the condensate 
fraction and the internal energy. In particular, because 
the energy E\ also varies as a function of T and N in our 
approach, we have to solve E\ and n with the conditions 
that fi cannot exceed E\ at any temperature. Further, 
since we take the parameter Na/aho ~ 1 which corre- 
sponds neither to the TF limit nor to the noninteracting 
limit for the semiclassical point of view, we cannot ex- 
pect jj, from the semiclassical approximation. Then, we 
find that the interaction effect shifts the transition tem- 
perature To given by /i = Eq at T = T c about 2 percent 
toward the lower temperature compared with the nonin- 
teracting model including the finite size effect. 



4 



ACKNOWLEDGMENTS 



I thank T. Fujita for careful reading. 



[1] M. H. Anderon, J. R. Ensher, M. R. Matthews, C.E. 

Wieman, and E. A. Cornell, Science 269, 198 (1995). 
[2] C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. 

Hulet, Phys. Rev. Lett. 75, 1687 (1995). 
[3] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van 

Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, 

Phys. Rev. Lett. 75, 3969 (1995). 
[4] M. -O. Mewes, M. R. Andrews, N. J. van Druten, D. M. 

Kurn, D. S. Durfee, and W. Ketterle, Phys. Rev. Lett. 

77, 416 (1996). 
[5] J. R. Ensher, D. S. Jin, M. R. Matthews, C. E. Wieman, 

and E. A. Cornell, Phys. Rev. Lett. 77, 4984 (1996). 
[6] W. Ketterle and N. J. van Druten, Phys. Rev. A 54, 656 

(1996). 

[7] H. Haugerud, T. Haugset, and F. Ravndal, Phys. Lett. 

A 225, 18 (1997). 
[8] E. P. Gross, Nuovo Cimento 20, 454 (1961), J. Math. 

Phys. 4, 195 (1963). 
[9] L. P. Pitaevskii, Sov. Phys. JETP 13, 451 (1961). 
[10] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 

Rev. Mod. Phys. 71, 463 (1999). 
[11] A. Griffin, Phys. Rev. B53, 9341 (1996). 
[12] S. Giorgini, L. P. Pitaevskii, and S. Stringari, J. Low. 
Temp. Phys. 109, 309 (1997). 



[13] D. S. Jin, M. R. Matthews, J. R. Ensher, C. E. Wieman, 

and E. A. Cornell, Phys. Rev. Lett. 78, 764 (1997). 
[14] S. Giorgini, L. P. Pitaevskii, and S. Stringari, Phys. Rev. 

A 54, 4633, (1996). 
[15] A. L. Fetter, Ann. Phys. (N.Y.), 70, 67 (1972). 
[16] L. D. Landau and E. M. Lifshitz, Statistical Physics, 3rd 

ed. (Pergamon Press, 1980). 
[17] G. Baym and C. Pethick, Phys. Rev. Lett. 76, 6 (1996). 
[18] M. Edwards and K. Burnett, Phys. Rev. A51, 1382 

(1995). 

FIGURE CAPTIONS 



FIG. 1. Chemical potential fi/ksTo as a function of T/T . 
The solid line refers to the noninteracting model. Solid circles 
show the calculated values of the interacting case. 

FIG. 2. Condensate fraction N/No as a function of T/Tq. 
The solid line is drawn with Eq. ( [l.l| ) and the dotted line is 
the thermodynamic limit. Solid circles are the same as in Fig. 



FIG. 3. The same as Fig. [2] just below the transition tem- 
perature. Open circles are numerical results of the noninter- 
acting case. 

FIG. 4. Internal energy E/Nk B T as a function of T/T . 
The solid line is the noninteracting case. Solid circles are the 
same as in Fig. |l| 
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